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Abstract 

We examine the phenomenological consequences of a supersymmetric bulk in the 
scenario of large extra dimensions. We assume supersymmetry is realized in the bulk 
and study the interactions of the resulting bulk gravitino Kaluza-Klein (KK) tower of 
states, with supersymmetry breaking on the brane inducing a light mass for the zero- 
mode gravitino. We derive the 4-d effective theory, including the couplings of the bulk 
gravitino KK states to fermions and their scalar superpartners. The virtual exchange 
of the gravitino KK states in selectron pair production in polarized e~^e~ collisions is 
then examined. We find that the leading order operator for this exchange is dimension 
six, in contrast to that of bulk graviton KK exchange which induces a dimension eight 
operator at lowest order. The resulting kinematic distributions for selectron production 
are dramatically altered from those in D = 4 supersymmetric scenarios, and can lead 
to a enormous sensitivity to the fundamental higher dimensional Planck scale, of order 
20 - 25 X 
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1 Introduction 



There has been much interest recently in the framework, proposed by Arkani-Hamed, Di- 
mopoulos, and Dvah (ADD) ||l], ^, which resolves the hierarchy problem by exploiting the 
geometry of spacetime. In this scenario, the fundamental scale of gravity in a higher D = 4+6 
dimensional spacetime is assumed to be of order the electroweak scale ~ 1 TeV. The appar- 
ent weakness of gravity in our 4-dimensional world originates from the large volume of the 
additional 6 spatial dimensions. The 4-dimensional Planck scale is no longer a fundamental 
scale, leaving the electroweak scale as the ultraviolet cut-off of the low-energy effective the- 
ory. The gauge hierarchy is thus effectively eliminated and reduced to the more tractable 
problem of stabilizing the higher dimensional radii In this scenario, gravity propagates 
throughout the higher dimensional volume, known as the bulk, whereas the Standard Model 
(SM) fields are confined to a 3-dimensional brane, or wall. 

In this theory Gauss' Law relates the Planck scale of the 4-dimensional theory, Mp, 
to the fundamental scale of gravity, M^, through the volume of the compactified dimensions 
Vs via 

Mj, = VsMl+', (1) 

where Mp = 1.2 x 10^^ GeV is the 4-d Planck scale. Setting Mp, ~ 1 TeV then determines 
the compactification radius Rc of the extra dimensions, with the exact relationship being 
set by the geometry of the compact dimensions. Assuming that the extra dimensions are 
of toroidal form and are all of equal size, we have Vs = {2t:RcY . then ranges from a 
sub-millimeter to a few fermi for 5 = 2 to 6. The case of 5 = 1 is excluded as it predicts 
corrections to Newtonian gravity at distances comparable to those in the solar system. A 
similar scenario can be realized in string theory where the string scale plays the role of 
the higher dimensional fundamental scale [Q, with the string scale acting as the ultraviolet 
cut-off of the theory. 
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Proposals for the localization of SM matter and gauge fields to a 3+1 dimensional wall 
have been made in the context of topological defects of higher dimensional field theories [^]. 
Such localization can occur naturally in string theory via D-branes where the SM particles 
are represented by open strings whose ends lie on the D-brane, while gravitons, which carry 
no gauge charges, may propagate in the bulk and correspond to closed strings 0, 0]. 

Since this scenario modifies gravity at the electroweak scale, it is natural to expect 
the emergence of new phenomena at the TeV scale which may reveal itself in experiments 
and lead to signatures very different from SM predictions. Upon compactification, the bulk 
graviton expands into a Kaluza-Klein (KK) tower of states, referred to as a bulk graviton KK 
tower, which interact with the SM fields on the brane. Collider signals for the graviton KK 
states have been studied by various authors P, |^, who have considered the virtual exchange 
of bulk graviton Kaluza-Klein towers, processes which radiate gravitons into the bulk, and 
stringy excitations of the Standard Model particles. Data from the Tevatron, LEPII, and 
HERA |jlO| presently constrain Md > 1 TeV for all values of 6, while the LHC and a future 
high energy e~^e~ Linear Collider are expected to probe fundamental scales in the 5 — 9 TeV 
range. Astrophysical and cosmological considerations place stringent bounds, of order 
Md ^ 100 TeV, for the case of 5 = 2; these limits weaken substantially to Md ^ 1 TeV 
for higher values of 6. Mechanical experiments have tested the inverse-square nature of the 
gravitational force law down to distances of 150/xm [|T2| for the case of 5 = 2. This scenario is 
thus consistent with all data as long as the number of extra dimensions is greater than one, 
with the case oi 6 = 2 being disfavored in terms of being relevant to the hierarchy problem. 

While the original motivation for the ADD scenario was to solve the hierarchy problem 
without the introduction of low-energy supersymmetry (or technicolor), one still might ask 
whether supersymmetry plays a role in such a scenario. Clearly, bulk supersymmetry is 
not in conflict with the basic assumptions of the model. In fact, various reasons exist for 



2 



believing in a supersymmetric bulk, not least of which is the motivation of string theory. As 
discussed above, D-branes of string theory provide a natural mechanism for the confinement 
of the SM fields. If string theory is the ultimate theory of nature then the proposal of 
ADD might be embedded within it with a supersymmetric bulk, i.e., a bulk supporting a 
supersymmetric gravitational action. In addition, extra dimensional geometries have been 



shown to provide novel methods for breaking supersymmetry ||T3|, |14|, the possibility that 
a supersymmetric bulk might provide a source for a tiny cosmological constant has been 
discussed ||15| , and supersymmetry might also serve as a mechanism for stabilizing the bulk 
radii. Supersymmetry has also been considered in the context of warped extra dimensions 
present in the Randall- Sundrum scenario of localized gravity |16| . 



In this paper, we investigate the consequences of a supersymmetric bulk in the ADD 
scenario. If bulk supersymmetry remains unbroken away from the brane, then it is natural 
to ask what happens to the superpartners of the bulk gravitons, the gravitinos. The bulk 
gravitinos must also expand into a Kaluza-Klein tower of states and induce experimental 
signatures. Up to now, the phenomenology of such gravitino KK states has been unex- 
plored. The existence of a light graviphoton, which is present if there is no orbifolding, in 
supersymmetric ADD has been shown to alter the non-supersymmetric ADD phenomenology 



Tq , and hence we also expect large effects from the gravitino sector. Here, for concrete- 
ness, we examine the virtual exchange of the bulk gravitino KK tower in superparticle pair 
production. 

We now outline the approach taken in this paper. We study the compactification of 
the gravitino sector of a generic supergravity theory living in the bulk and work out the form 
of the effective action describing the free part of the reduced action on the 3 + 1 dimensional 
brane to which the fields of the Minimal Supersymmetric Standard Model (MSSM) are 
assumed to be confined. We assume that supersymmetry on the wall is broken, giving rise 
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to masses for the zero-modes of the Kaluza-Klein tower resulting from the reduction of 
the higher dimensional gravitinos, and shifting the masses of all higher modes. We make no 
special assumptions about the nature or origin of confinement of the MSSM fields to the wall. 
We then derive the couplings of the Kaluza-Klein modes of the higher dimensional gravitinos 
to fermions and their scalar superpartners on the wall, yielding an effective 4-d theory which 
we then use to study the corrections to collider signatures of certain processes in the MSSM. 
For purposes of illustration, we work with a 10 dimensional theory; the generality of our 
results should be apparent (up to considerations of representations of fermions in various 
dimensions) . 

The gravitino zero-mode acquires a mass from the spontaneous breaking of super- 
symmetry on the brane, with this mass being proportional to ~ A|jjgy/Mp, where Asusy 
represents the supersymmetry breaking scale in a generic model, noting that the zero-mode 
gravitino couples with the usual Mp^ strength. This familiar 4-d expression can also be seen 
to arise from the volume factor after integrating over the 6 extra dimensions and using Eq. 
(1). Various SUSY breaking scenarios thus yield different predictions ranging from ultra- 
light gravitinos from weak-scale supersymmetry breaking models, light gravitinos in gauge 
mediated SUSY breaking, and heavier gravitinos in models of gravity and anomaly mediated 
SUSY breaking. The possibility of light gravitinos is an interesting one and their generic col- 



lider [0 and astrophysical ||19[ implications have been studied, with direct collider searches 
yielding a bound of ^ 10~^ eV [^. In particular, the case of gauge mediated SUSY 
breaking (GMSB) has been intensively studied |2T[], as a light gravitino modifies the standard 
collider search techniques for supersymmetry and provides a good candidate for warm dark 
matter p2|. In this paper, we work in the context of gauge mediated supersymmetry as it 
naturally affords a light gravitino state, however our results apply to any supersymmetric 
model with a light gravitino, including weak-scale breaking models. 
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We focus on the effects of the virtual exchange of the bulk gravitino and graviton KK 
tower states in the process e"'"e~ e+e^ at a high energy Linear Collider (LC). This process 



is well-known as a benchmark for collider supersymmetry studies , as the use of incoming 
polarized beams enables one to disentangle the neutralino sector and determine the degree 
of mixing between the various pure gaugino states. The effects of the virtual exchange of a 



bulk graviton KK tower in selectron pair production has been examined in pj] for the case 
of non-supersymmetric large extra dimensions. Here, we will see that the introduction of 
bulk gravitino KK exchange greatly alters the phenomenology of this process by modifying 
the angular distributions and by substantially increasing the magnitude of the cross section. 
We find that the leading order behavior for this process is given by a dimension-6 operator, 
in contrast to the dimension-8 operator corresponding to graviton KK exchange. This yields 
a tremendous sensitivity to the existence of a supersymmetric bulk, resulting in a search 
reach for the ultraviolet cut-off of the theory of order 20 — 25 x ^/s. 

Our paper is organized as follows: In Section 2 we present the Kaluza-Klein decom- 
position of the gravitinos in such a model. In Section 3 we use the general Noether technique 
to couple the gravitino Kaluza-Klein states to matter fields, and in Section 4 we present the 
phenomenology of this scenario, pointing out some generic features of such models. Our 
conclusions are given in Section 5. Various details are relegated to the appendices. 

Our notation in this paper is as follows: We assume space-time possesses 3 + 1 + 
6 dimensions, with matter confined to a 3 + 1 dimensional brane, and pure supergravity 
propagating also in the extra S dimensions. Letters from the Greek alphabet are used to 
denote curved (world) indices while those from the Latin alphabet denote flat (tangent space) 
indices. Hatted symbols range over the full 4 + 6 dimensions; barred ones are restricted to 
the 6 extra dimensions, while those without any decoration live on the 3-1-1 dimensional 
brane. Coordinates over the entire = 4 + 6 dimensional manifold split as z^^ = {x^, y^). 
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Our Minkowski metric convention is 77-^ = (+1, —1, . . . , —1). 

2 Kaluza-Klein Excitations of the Gravitino 

We begin by reminding the reader of some general observations about string theory and su- 
pergravity in ten and eleven dimensions. Recall that type IIA (ten dimensional) supergravity 
with two supersymmetries of opposite chirahty can be deduced via dimensional reduction 
of the unique 11 dimensional (N=l) supergravity theory. The type IIA superstring theory 
whose low energy effective action gives rise to the type IIA supergravity admits Dp-branes 
with p even, and is S-dual to eleven dimensional M-theory |^ , with the eleven dimensional 
supergravity as its low energy limit. M-theory is believed to unify the five known string 
theories. The type IIB string theory (which admits Dp-branes with p odd) is T-dual to the 
type IIA theory. 

D-branes are extended string theoretic objects on which open strings terminate [0]. 
Only closed strings can propagate far away from the D-brane, which on a ten dimensional 
background is described locally by a type II string theory whose spectrum contains two 
gravitinos (their vertex operators carry one vector and one spinor index). The D-brane 
introduces open string boundary conditions which are invariant under just one supersymme- 
try, so only a linear combination of the two original supersymmetries survives for the open 
strings. Since open and closed strings couple to each other, the D-brane breaks the original 
N = 2 supersymmetry down to = 1. The low energy effective theory will then be a 
D = 10, N = 1 supergravity theory with a single Majorana-Weyl gravitino which couples to 
a conserved space-time supercurrent. For our model, we assume that the Standard Model 
fields are confined to a 3-brane, with pure supergravity living in the bulk of space-time. 
Attempts at constructing a Standard Model on D-branes can be found in pB| . 
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Pure supergravity in four dimensions contains as physical fields the vierbein e and 



the gravitino \1'^ p7| , P8| . In the free field limit, the gauge action of supergravity reduces to 
the standard Fierz-Pauli action (the linearized part of the Einstein-Hilbert action), together 
with the Rarita-Schwinger action, which are the unique ghost free actions for a spin 2 and 
spin I field [^, 0. 

The spin of a field ^{z) is given by a representation of the tangent space group at a 
fixed point z Under the above compactification, the metric (^^^ decomposes in four 

dimensions as a metric (7^1,, while the components Qf^p with one index ranging over the extra 
dimensions are seen in four dimensions as a vector, and g/ip transforms as a scalar. To study 
the behavior of spinors on a general manifold, we must introduce a set of orthonormal basis 
vectors (the vielbein) e^^{z) in the tangent space of the manifold at the point 2, with fi 
transforming as a vector index under general coordinate changes, while m transforms as a 
scalar under such diffeomorphisms, but as a local Lorentz vector index under SO{D — 1,1), 
i.e., the Lorentz group acting in the tangent space at the point z, satisfying |3^, |3^ 



r^^- = /'^e^e/ , g/^o = r/^^e>e> . (2) 

Associated with the local Lorentz symmetry is a gauge field for local Lorentz transforma- 
tions, the spin connection Ufimn, transforming as a covariant vector under general coordinate 
transformations. 

We now assume the vacuum of space-time to be of the form x T^, where 
is four dimensional Minkowski space-time and = x . . . x , the direct product of 
six dimensions each compactified on a circle. This product form of the vacuum insures 
four dimensional Poincare invariance. Physical fluctuations of space-time can be treated in 
perturbation theory by expanding the metric around this vacuum. 

Consistent with the symmetry of the vacuum state, we assume the vielbein takes the 



form 




e/(t/) 



(3) 



where we have used the local SO (9,1) gauge symmetry to set the E'^'* components to zero. 
Here ej^ is the vielbein on the compactified space, are the massless gauge fields of the 
symmetry group of the internal space (here t/(l)®), and the are the Killing vectors 
associated with this symmetry. 

The gravitino kinetic term of the Lagrangian is 



with E being the determinant of the vielbein in ten dimensions, V^^^ the antisymmetric 
product of three F matrices defined by 



and '^fi a Majorana-Weyl vector-spinor. We have suppressed the spinor indices for notational 
clarity. 

We use F'^ to denote Dirac gamma matrices for the full D dimensional space, and 
to denote them on the 3 + 1 dimensional brane. As usual, the Dirac algebra (in the tangent 

space) is spanned by a set of constant matrices (which in 10 dimensions are 32 x 32) 



(4) 




(5) 



mn 



(6) 




(7) 



satisfying 



{F^(z),F^(z)} = 2gf''{z). 



(8) 
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A convenient representation of the algebra in ten dimensions which makes manifest 
the four dimensional decomposition is given in Appendix ^. The S0(9,l) generators are 

j.mh ^ ^ [r'^^r^] . (9) 

The first four matrices in the set ( p3D furnish a reducible representation of S0{3, 1). Taking 

account of the identity acting in the internal 2^ = 8 dimensional spinor space, we see that 
a -D = 10 Dirac spinor decomposes into 8 D = A Dirac spinors. A D = 10 Majorana spinor 
would decompose into 8 Majorana spinors in four dimensions, and a D = 10 Majorana- Weyl 
spinor would decompose into four Majorana spinors, since the D = 10 chirality condition 
rii\l/ = ±\E' (where Tu = iTq . . . Fg plays the role in ten dimensions of 75 in four dimensions) 
would pair-wise relate half the degrees of freedom, leaving four D = 4 Majorana spinors. 
There will also be 24 Majorana spin-| fields. After this reduction we recover the standard 
definition of the four dimensional spinor generator. The metric in ten dimensions decomposes 
into one spin-2 graviton in four dimensions, 6 vectors, and 21 scalars. 

In ten dimensions, a Majorana- Weyl spinor contains 16 real components. To see 
this, note that in ten dimensions, the Dirac matrices are 32 x 32, so a Dirac spinor would 
contain 32 complex components. The Majorana condition reduces this in half by imposing 
a reality condition. Finally, the Weyl condition reduces this again, for a final total of 16 
real components. The equations of motion for this spinor imply that not all of the remain- 
ing components are independent, reducing further the independent propagating degrees of 
freedom. 

For simplicity we assume the compactification radii of all six dimensions to be the 
same, though the generalization is obvious. Physically this compactification amounts to the 
identification 

= + 2nnR, , (10) 
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with n an arbitrary integer and Rc the common radius of the compactified dimensions. The 
condition on the fields then becomes, recaUing our notation z^^ = (a;^, yP"), 

/ + 27rnR,) = ^^(x^ /) . (11) 

The transformation properties of the Rarita-Schwinger field \E'^ are the product of a 
vector and a spinor. The covariant derivative is given by 

^p. = dfi + —LOfirrm'^"^"'- (12) 

where S™" are the SO (9,1) generators given in (||), and uixmh is a sum of the standard spin 
connection and terms quadratic in gravitino fields with supercovariantization giving rise 
to four fermion terms. 

Using the vielbein to translate from curved to fiat indices, we can rewrite the action 
d^), after linearizing the spin connection, in the form 

E-'£. = ^ r'^"^ 9^v|>, , (13) 

where 

^^ = i?/v^., (14a) 

dn = Ef^'dc (14b) 

The linearized action possesses the local Abelian gauge symmetry 

S-^rh = d^e, (15) 

with e a Majorana-Weyl spinor. This is the same linearized gauge invariance one would 
expect of the fermionic part of the gauge action of supergravity. The above invariance 
follows from the realness of £ and the total antisymmetry of r"^"-P, 



10 



To make the Kaluza-Klein decomposition we expand the Rarita-Schwinger field 
in eigenf unctions of the compactified space (of dimension 6, which for now is taken to be 
arbitrary), assumed here to be the torus T^, with volume Vs = (2tiRcY . The expansion is 



The four dimensional fields are seen to arise as the coefficients in this expansion. 

We take as our starting point the linearized action ([T3|). We now substitute ([T6| ) to- 
gether with (^) in the action (|13|), then decompose the 10 dimensional indices as m = (m, m), 
with the indices m and fh transforming as vectors under 5*0(3, 1) and S0{5) respectively. 
The Rarita-Schwinger field then splits as 

^rh = (^m,^m), (17) 

where 

^m = EJ^^, (18a) 
<frn = Ej<i/i, = e/^^. (18b) 
Introducing the shorthand notation 

we find 

oo oo 



s=-oo t=-oo 



+ $?,r'^"^'(9„<i>y + $?,r™p(9„<i>J)) e"+^ 
+ (^¥j:^''p{dne^)¥^ + $lr'""p"(a^e^)$J 

+ l>^r'^"P(a^,e^)<l>; + ¥j:'^''P{dne'')%) e 
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(20) 



with e being the determinant of the four dimensional vierbein. The first set of terms propor- 
tional to e""*"^ and with the derivative acting on the fields and $^ are the kinetic terms 
since the derivative is taken with respect to coordinates in M*^, while the remaining terms 
proportional to e° with the derivative acting on the exponential will give rise to the mass 
terms of the spin | and | fields in the effective four dimensional action. It is important 
to note that fields with indices ranging over the 3 + 1 dimensions ($m) appear in this La- 
grangian coupled to those with indices in the extra 6 dimensions in both the kinetic 

and mass terms. This mixing prevents us from interpreting $m as the four dimensional 



Rarita-Schwinger field. We introduce the field redefinition |36, 37 



xiix) ^ + ^r„r$^:(x) , (21) 

which will be identified as the four dimensional Rarita-Schwinger field. In what follows, 
we shall systematically drop all terms coupling only fields with vector indices in the extra 
dimensions (these would give rise to spin | terms in the four dimensional effective theory), 
and concentrate on the gravitino (spin |) part of the effective theory. The field redefinition 
([2T| ) completely decouples the kinetic term into spin | and spin | pieces, but the mass term 
still contains couplings between spin | and | fields. The mixing with spin ^ in the mass 
term is due to the fact that the field Xm transforms as the 



1 1 
2' 2 



io)e(o.l 



i,i)e(i.i)®(io)®(o.i) (22) 



representation of S0(l,3). To project out the (|,0) © (0, |) part we impose the S0(3,l) 
invariant gauge condition 

r-x^(x) = 0, (23) 

leaving only a spin | field. 
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To evaluate the derivatives in the mass terms of (^), we note that they are taken 
over coordinates in the extra dimensions, and since y ranges over the same coordinates, they 
yield 

dne^ = -^Ue^ (24) 



Putting together (pol), (piD and (P4|), we see that the spin | part of the Lagrangian takes the 
form 



where xtn — (Xm)W ® Is? the S™" are given by restricting the indices of the S0(9,l) 
generators (y) to the first 3 + 1 dimensions, and we have summed over n ranging over the 
extra dimensions. 

The only dependence on the coordinates of the extra dimensions in (^) is isolated 

in the exponential e°'~^^. Integrating over the extra dimensions gives the four dimensional 
effective Lagrangian 

£e//(x) = / dyi...dysJ^{x,y). (26) 



The integrals in ( PB| ) can be performed using the orthogonality of the exponentials, with the 
result 

/ dyi... dys exp - — (s - t) ■ y ='^5 17 , (27) 

where = {2ttRc)^ is the volume of the compactified extra dimensions. The effective 
Lagrangian defined by (|26|) can be written as a sum over Kaluza-Klein states 

oo 

^e//(x) = (28) 
s=—oo 
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with the Lagrangian for each level 

e-'Liff{x) = lxi{x)r--ndnX;{x)) - zxl(x)S™n;^r%)x;», (29) 

and the fields Xm describing an individual Kaluza-Klein mode. In (|29|), the r™"^ and S™^ 
are still 32 x 32. We note that masses of the Kaluza-Klein excitations of the gravitino are the 
same as for the gravitons. To arrive at a proper gravitino mass term, we must diagonalize 
the term T"'Sn in the internal space. First we note that in the representation ( |53D 

Ymnp ^ ^mnp ^ ^ 

(30) 

= cr ® Is , 

with the 7™"'P and a'"" being products of 4 x 4 four dimensional Dirac matrices. We decom- 
pose the fermions into four dimensional ones Xa ^a,j with a = 1 ... 32 and a = 1 ... 4 
spinor indices in ten and four dimensions respectively, and j = 1 . . .8 serving as an internal 
index, and then use xfn — (Xm)W^®l8- Making a unitary transformation on the fields allows 
us to diagonalize T^Sn without affecting the kinetic term since in the chosen representation, 
it is already diagonal in the internal space. After diagonalization in the internal subspace, 
the Lagrangian can be written as 

e-' L%f{x) = [l^'n.\x)r'''idn4''i^)) + oo'^\x)a-^^'m^;^^{x)^ , (31) 

with the mass eigenvalue coming from the diagonalization and is given by mi = {—ly^^^. 

The fields associated with the negative mass eigenvalues can be redefined to remove this 
sign, however, care must be taken with the Feynman rule for the coupling of the gravitino 
to matter. The 7^ in the mass term can be removed via a chiral rotation of the fields, which 
introduces an extra factor of i into the mass without affecting the kinetic term. The sum 
in (|31| ) runs over the four Majorana vector-spinors in four dimensions. We have applied 
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the Majorana-Weyl condition in ten dimensions, which, as previously discussed, yields four 
Majorana spinors after the decomposition into four dimensions. Generally, the masses of 
the four gravitinos at each Kaluza-Klein level can be shifted by supersymmetry breaking 
effects on the brane. When we consider the phenomenology of such models, we will assume 
that the = 4 supersymmetry is broken at scales near the fundamental scale with 
only = 1 supersymmetry surviving down to the electroweak scale. The phenomenological 
contributions from the heavy gravitinos associated with the breaking of the extended super- 
symmetry near the fundamental scale will be highly suppressed, due to the large mass for 
these individual excitations. 

Using the identity 

the compactified action for each gravitino can be put in the standard form for the action of 
a massive spin | particle in four dimensions 



The mass term appearing in ( p3| ) has the same form as the mass generated by spontaneous 
breaking of supersymmetry. The massless limit of the propagator associated with this La- 
grangian contains singular terms, but these terms are all proportional to or g^, and since 
ujm couples to a conserved current (the Noether current associated with supersymmetry, 
hence the interpretation as the gauge field of supersymmetry), their contributions vanish. 
We assume that supersymmetry on the brane is broken, giving a mass to the lowest lying 
Kaluza-Klein states (s = 0) in (0), and additively shifting the masses of the higher modes 
by rrifi — >■ rrifi + mQ. A natural mechanism for breaking supersymmetry which ensures a light 
gravitino is gauge- mediation |]2T]]. In these models, the gravitino is the lightest supersym- 



metric particle, and if R parity is a good symmetry, it is the terminus for the the decay 
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chains of all superparticles. The lightness of the gravitino simplifies the summation over the 
Kaluza-Klein states (see Appendix C) leading to results which are essentially independent 
of the mass of the zero-mode. 

The Lagrangian (^3|) implies the field equation 

e"^"^^57n5.<(x) + -f [7™,7lO) = (34) 

obeyed by each Kaluza-Klein excitation. Contraction of ( p4|) with 7^ yields 



r^oo^ = 0, (35) 

since by taking the divergence of (|34D , we have [j^, 7"] cj*'-' = 0. This shows that the field 
equation (0) indeed describes a particle of spin | with no spin ^ admixture ||3^. Each 
vector component (m) of the Rarita-Schwinger field ui^ satisfies a Dirac equation. 

We note that the kinetic part of the compactified action still possesses the local gauge 
symmetry arising from the decomposition of (|T3p, but the mass term breaks this symmetry, 
thus allowing us to invert the free quadratic operator in the Lagrangian to find an appropriate 
propagator, which we derive in Appendix B. 

3 Gravitino Couplings to Matter 

In this section we discuss the coupling of fermions and scalars to gravitinos. We begin by 
formulating a globally supersymmetric theory, and then proceed to gauge this supersymme- 
try via the Noether procedure, deriving a locally supersymmetric Lagrangian yielding the 
coupling of the fermions and scalars to the graviton and gravitino. The Noether procedure 
provides a systematic technique for deriving an action with a local symmetry from an action 
possessing a global one PSI, 153]. Gauging the rigid supersymmetry transformations yields the 
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coupling of the Rarita-Schwinger field, which is the gauge field of the local supersymmetry, 
to the matter multiplet. 

The 7 matrices are, in the chiral representation, given by 

[ . 

7^^ = , (36) 

yaf" 

with 

a^ = (l,a), a^ = {l,-a), (37) 

and a being the Pauli matrices. A Dirac spinor is constructed as follows: 



where we have introduced dotted and undotted indices and ipL/R are, in the massless 



limit, states of definite helicity each satisfying a Weyl equation. 

We begin with the following off-shell Lagrangian in 3 + 1 dimensions (two copies of 
N = 1, D = 4, chiral supermultiplet) 

£ = (5^$t)(5^^$) + i^7^9^^ + 3^3^, (39) 

which describes a free massless fermion, a multiplet of complex scalar fields, and a set of 
auxiliary fields, with the complex scalar fields written as 

Here, F and G are non-dynamical (non-propagating) in the sense that their equations of 
motion are algebraic constraints that allow us to eliminate them from the action. They have 
been introduced to ensure that the number of bosonic and fermionic degrees of freedom match 
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off-shell and that the algebra closes without use of the equations of motion. After elimination 
of F and G, the supersymmetry algebra no longer closes and demonstration of invariance 
requires the use of the equations of motion for the remaining fields. The supersymmetry 
transformation parameter B forms a Dirac spinor built from left- and right-handed Weyl 
parameters 



Decomposing the Lagrangian (^) into a sum of left- and right-handed parts yields 



which is a sum of two Wess-Zumino actions, one each for left- and right-handed multiplets. 

The Lagrangian (^) is invariant, up to a total divergence which does not contribute 
to the action (at least for topologically trivial field configurations), under the following set of 
supersymmetry transformations, which we write in Weyl component form, forming a closed 
algebra on the fields 




(41) 



lileft = (9^0*)(5» + ^V^a(a^)""9^V^a + F*F 



bright = {d,p*){d^p) + IX'^i^ncad.r + G*G 



(42a) 




'a 




.a 




(43) 



5,F 



6cG 
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We now take the fields in the left- and right-handed multiplets to transform under 
the same supersymmetry transformations (yielding = 1 supersymmetry), with O a Ma- 
jorana spinor, and apply the Noether procedure to derive the coupling of the supergravity 
multiplet to the matter fields. The gravitino will appear as the gauge field of local super- 
symmetry. Supersymmetry breaking will manifest itself through the appearance of a mass 
for the gravitino. 

To gauge the global symmetry, we let the supersymmetry transformation parameter 
become space-time dependent O — > Q{x). The Lagrangian will no longer be invariant, 
but will vary as 

= J^d^Q, (44) 

with J'^ the global supersymmetry current. We must now add terms to the Lagrangian and 
the supersymmetry transformation rules until we restore the invariance of the Lagrangian 
(up to a total derivative). We add a term coupling the local supersymmetry gauge field to 
the symmetry current 

£1 = cfi^ J^ (45) 

with a Majorana vector-spinor, which we expect to become the gravitino, and hence must 
transform as 

= -d^Q (46) 

hi 

to leading order in n in the supergravity theory (here k = \^8itGn). The local variation 
vanishes to 0(k°) if c = — I, but not to order k. Iterating this process to higher orders in k 
and covariantizing with respect to gravity, we arrive at a locally supersymmetric Lagrangian 
at order k^. The result consists of the Einstein-Hilbert and the Rarita-Schwinger actions 
together with the original action (p9D minimally coupled to gravity, plus a term coupling 

19 



spinors, scalars and gravitinos, and higher order four point terms. The term couphng scalars, 
spinors and gravitinos, and minimally coupled to gravity, is 

Lj = -^\e\{{d,<^L)Q,'y^^Y^L + {d,<l>R) n^YY^n] + h.c, (47) 

with Qi, a Majorana vector-spinor. Expanding |e| to leading order in the vierbein yields the 
Feynman rule displayed in Fig. |l|, with the obvious generalization to the hermitian conjugate 
piece. 




Figure 1: Feynman rule for the gravitino, fermion, scalar coupling. Pl,r represent the 
standard projection operators. 



4 Phenomenological Analysis and Numerical Results 

We are now ready to apply our results to a phenomenological analysis. For purposes of 
illustration we focus on selectron pair production in high energy polarized e~^e~ collisions. 
As discussed in the introduction, this process provides a valuable tool within the MSSM for 
determining the composition of the mixed neutralino states, iii terms of the various pure 

SU(2) and U(l) wino and bino components, W^, . Here, we investigate how the existence 
of supersymmetric extra dimensions modifies this reaction. 
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^R^R 


e^e+ 


s-channel , Z ,Gn 
t-channel W ,B 




t-channel B , 


S-channel 7,2', 




s-channel 7 , Z , 


t-channel B , 




s-channel 7,2', 
t-channel B , 



Table 1: The diagrammatic contributions to individual scattering processes for polarized 
electron beams. A blank box indicates that there are no contributions for that polarization 
configuration. 



The tree level processes contributing to selectron pair production in the presence 
of a supersymmetric bulk are presented in Fig. |^. In addition to the standard 7, Z s- 

channel exchange and B^, t-channel contributions present in the MSSM, we now have 
contributions arising from the s-channel exchange of the bulk graviton KK tower and the 
t-channel exchange of the bulk gravitino KK tower. There are no u-channel contributions 
due to the non-identical final states. The contributions from neutral higgsino states are 
negligible due to the smallness of the Yukawa coupling. The diagrammatic contributions 
to the individual scattering processes for left- and right-handed selectron production with 

initial polarized electron beams are summarized in Table |l|. Note that the exchange 

only contributes to the process e^e"*" —* e^e^, and that the t-channel gravitino and the 5° 

contributions are isolated in the reaction e^^e+ '^lr'^rl- 



Our amplitudes for the standard MSSM contributions to the reaction e~ {p2)e~^{pi) — >■ 
e~ {k2)e~^{ki) (with the direction of the charge flow as indicated in Fig. ^ reproduce the re- 
sults in The unpolarized matrix element for the case of massive gravitino KK exchanges 
is 

^ = T E T±\^iPih,l-^P''''^%lre{p2) , (48) 

n " 
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Figure 2: Processes contributing to scalar electron pair production in polarized e^e col- 
lisions. The state h^^ represents the bulk graviton KK tower, the bulk gravitino KK 

tower, and 5°, correspond to the unmixed MSSM gaugino states. 
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where the sum extends over the gravitino KK modes and we recall that k = y/SnGN = Mp^ 

is the reduced Planck scale. P"'^"^ represents the numerator of the propagator for a Rarita- 
Schwinger field of mass and is given in Appendix ^. The mass splitting between the 
evenly spaced bulk gravitino KK excitations is given by 1/Rc, which lies in the range 10~^ 
eV to few MeV for 5 = 2 to 6 assuming Md ~ 1 TeV; their number density is thus large at 
collider energies. The sum over the KK states can then be approximated by an integral which 
is log divergent for 5 = 2 and power divergent for 6 > 2. We employ a cut-off to regulate these 
ultraviolet divergences, with the cut-off being set to Ac, which in general is different from 
Mjy, to account for the uncertainties from the unknown ultraviolet physics. This approach 
is the most model independent and is that generally used in the case of virtual graviton 
exchange f^]. In practice, the integral over the gravitino KK states is more complicated than 
that in the case with spin-2 gravitons due to the dependence of the gravitino propagator on 
m^j. We find that the leading order term for \/\t\ ^ Ac results in the replacement (in the 
case of 5 = 6) 

in the matrix element; the structure of the summed gravitino propagator is thus altered from 
that of a single massive state. Hence the leading order behavior for gravitino KK exchange 
results in a dimension-6 operator! This is in stark contrast to graviton KK exchange which 
yields a dimension-8 operator at leading order. We thus expect an increased sensitivity to 
the fundamental scale Md in the case of a supersymmetric bulk. Our derivation of the 
summation over the bulk gravitino KK states in the matrix element is detailed in Appendix 



C. Lastly, the bulk graviton KK tower contributes an additional amplitude of the form 

2 



^ = T7(^ - - k2)AP2h''e{p,) . (50) 



In order to perform a numerical analysis of this process, we need to specify a con- 
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Crete supersymmetric model. We choose that of Gauge Mediated Supersymmetry Breaking 
(GSMB) as it naturally contains a light zero-mode gravitino. We specify a sample set of 
input parameters at the messenger scale, where the supersymmetry breaking is mediated via 
the messenger sector, and use the Renormalization Group Equations (RGE) to obtain the 
low-energy sparticle spectrum. We choose two sets of sample input parameters describing 
the messenger sector which are consistent with our model. The RGE evolution of these 
parameter sets is performed in ISAJET v7.51 EDI and results in the sparticle spectrum 



Set I : mg^ = 217.0 GeV, mg^ = 108.0 GeV, x° = (76.5, 141.5, 337.0, 367.0) GeV , 
Set II : mg^ = 210.5 GeV, mg^ = 104.5 GeV, x° = (110.5, 209.6, 322.5, 324.0) GeV , 

where Xi with i = 1,4 corresponds to the four mixed neutralino states. The first set of 
parameters yields a bino-like state for the lightest neutralino, whereas the second set results 
in a Higgsino-like state for Xi- These two sets are chosen so we can investigate the dependence 
of the kinematic distributions on the composition of the lightest neutralino state. Note that 
the ei and cr masses are essentially equivalent between the two sets and hence will not 
induce any kinematical differences in the distributions. In addition, these input parameters 
were selected in order to obtain a sparticle spectrum which is kinematically accessible to the 
Linear Collider; our results are essentially insensitive to the exact details of the spectrum. 
We stress that our analysis is purely phenomenological and that our conclusions do not 
depend on the physics inherent to GMSB. Except where noted, we perform our numerical 
analysis for the case of 5 = 6, following our above discussion of supergravity models. From 
here on, we refer to these two spectra as our D = 4 supersymmetric models. 

It is instructive to first examine the effects of each class of contributions to selectron 
pair production. This is displayed in Fig |^, which shows the angular distribution for the 

process e^e"*" — >• e^e^ with y/s = 500 GeV assuming 100% polarization of the electron 

beam; we show this particular reaction merely for purposes of demonstration. The bottom 
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curve represents the full contributions (s- and t-channel) from the 4-dimensional standard 
gauge-mediated supersymmetric model discussed above in the case where the Xi is bino-like, 
corresponding to parameter set I. Our numerical results for the MSSM case agree with those 
in the literature The middle curve displays the effects of adding only the s-channel 

contributions of the bulk graviton KK tower in the scenario of a non-supersymmetric bulk 
with Ac = 1.5 TeV. We see that there is little difference in the distribution between the 
D = A supersymmetric case and with the addition of the graviton KK tower, in either shape 
or magnitude. It would hence be difficult to disentangle the effects of graviton exchange from 
an accurate measurement of the underlying supersymmetric parameters using this process 
alone. The top curve corresponds to the full set of contributions from a supersymmetric 
bulk, i.e., our standard supersymmetric model plus KK graviton and KK gravitino tower 
exchange for the case of six extra dimensions with Ac = 1.5 TeV. Here we see that the 
exchange of bulk gravitino KK states yields a large enhancement in the cross section and 
a substantial shift in the shape of the angular distribution, particularly at forward angles, 
even for Ac = 3^/s. This provides a dramatic signal for a supersymmetric bulk! 

In Fie. H we explore the modifications to the angular distribution for Cp^e^ — > ^~r^r 
from a supersymmetric bulk as the value of 6 changes, using our parameter set I for demon- 
stration. In principle there are two competing effects which may modify the distributions: (i) 
the number of degenerate gravitinos in each KK level as a result of the reduction of fermions, 
versus (ii) the volume factor that appears in the density of states in the integral over the 
gravitino propagator which sums over the states in the KK tower, and its dependence on 
Ac. As discussed above, the number of degenerate gravitinos is reduced as the number of 
extra dimensions decreases; this results in a reduction of the cross section in the general case 
of extended low-energy supersymmetries. However, this effect does not modify our analysis 
since we have assumed the = 1 at low-energy. The second effect arises from the increase 
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in i?c, with Ac being held fixed, for smaller values of S as can be seen in Appendix C. This 
volume factor arises in the integral over the propagators for the bulk gravitino KK tower and 
is discussed in Appendix C. In this figure, we compare the event rates for this process for the 
cases 6 = 2 and 6, corresponding to the top and bottom curves, respectively. We note that 
the shape of the distribution differs in the two cases due to the form of the sub-leading terms 
in the integral over the propagator of the KK states. For the remainder of our analysis, we 
will display results only for the more conservative case of 5 = 6. 
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Figure 3: The angular distribution for e^e^ (^r'^r from the D = A supersymmetric model 
I, plus the addition of bulk graviton KK tower exchange, and with bulk gravitino KK tower 
exchange, corresponding to the bottom, middle, and top curves, respectively. 



Let us now study the variations in the distributions between the two different com- 
positions of the lightest neutralino. Figure ^ shows the angular distributions with 100% 
electron beam polarization for each helicity configuration listed in Table [l| for the two sets of 
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Figure 4: A comparison of the event rates for the cases with 6 = 2 and 6, corresponding to 
the top and bottom curves, respectively. 

parameters discussed above, with and without the contributions from supersymmetric extra 
dimensions. In each case, the sohd curve corresponds to the bino-hke case and the dashed 
curve represents the Higgsino-hke scenario. The top set of curves are those for a super- 
symmetric bulk with Ac = 1.5 TeV, while the bottom set corresponds to our two D = 4 
supersymmetric models, i.e., without the graviton and gravitino KK contributions. We note 
that the = 4 results agree with those in the literature . We see from the figure that in 
the process where the gravitino contributions are dominant, e^^e"*" — > e^e^, there is little 
difference in the shape or magnitude between the two Xi compositions. The use of selectron 
pair production in polarized e~^e~ collisions as a means of determining the composition of 
the lightest neutralino is thus made more difficult in the scenario with supersymmetric large 
extra dimensions. In what follows, we present results only for the bino-like Xi as a sample 
case; our conclusions will not be dependent on the assumptions of the composition of the 
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lightest neutralino. 

We now examine in Fig. ^ the total cross section as a function of center-of-mass 
energy for each helicity configuration. In each case, the bottom curve represents the D = 4 
bino-like supersymmetric model, while the remaining curves, from top to bottom, are for a 
supersymmetric bulk with Ac = 1.5, 3.0, 6.0 TeV. In some reactions, the results for Ac = 6.0 
TeV are indistinguishable from the D = 4 case. Here we can see the effects of unitarity 
violation as y/s approaches the value of the cut-off scale. Clearly, the new, as of yet unknown, 
ultra-violet physics will set-in at this point to regularize the cross section. 

Next, we present in Fig. |^ and || the number of events for the binned angular dis- 
tribution for each helicity configuration with 80% polarized electron beams for ^/s = 500 

GeV and 500fb^^ of integrated luminosity. In each case, the solid histogram corresponds to 
the D = 4 bino-like supersymmetric model, while the "data" points represent the addition 
of the bulk graviton and gravitino KK tower exchange for Ac = 1.5,3.0,6.0 TeV from top 
to bottom. As before, the contributions with Ac = 6.0 TeV are only distinguishable from 
the D = 4 results in the case of e^^e^ ^l^/j- The error bars on the "data" points are 
statistical only. We see that in most of the reactions, the case of A = 3.0 TeV leads to only 
a slight increase in event rate in each bin, whereas for e^^e"'" — e^e^, the t-channel bulk 

gravitino KK exchange is significant, leading to observable deviations from the D = 4 case 
even for Ac = 6.0 TeV. 

An interesting polarization asymmetry can be defined for the case of e^e^ and e^e^ 



28 



J3 
03 

r/l 


■0 
\ 
b 








b 
■0 



500 



100 
50 





1 ' ' ' ' 1 ' ' ' ' - 


A, = 1.5 TeV^^^f^T 




^ - ' 


(a) ei,e ^ e[, ^ 

Vs = 500 GeV - 
P^- = 100% 




(5 = 6 

1 , , , , 1 , , , , 



-1.0 



-0.5 



0.0 
COS £ 



0.5 



1.0 




-0.5 0.0 
COS B 



600 



■0 

\ 400 
b 



2000 



1000 




30 



S5 



to 

m 


V 

\ 

b 

T3 



' ' ' I ' ' 

(d) eae"" ^ Cl 



A„ = 1.5 TeV 




(e) e^e -» Cr 






^ y^ 

y yy''^ 
y yy 
y / 
y / 


/■ 
// 


y \, = 1.5 TeV 


// 
// 

// 










\\ 




— y'^ ^ 

-/ , , , 1 , , , , 1 





-1.0 



-0.5 



0.0 
COS S 



Figure 5: Angular distributions for each helicity configuration with supersymmetric bulk con- 
tributions for Ac = 1.5 TeV (top curves), and for the D = 4 supersymmetric models (bottom 
curves). The solid (dashed) curves correspond to a bino-like (Higgsino-like) composition of 
the lightest neutralino. 
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Figure 6: Total cross sections as a function of center of mass energy. The bottom curve 
corresponds to the D = 4 supersymmetric model, and the other curves are for a. 6 = Q 
super symmetric bulk with the cut-off as labeled. In figures (c) and (d), the Ac = 3 TeV 
curve is barely distinguishable from that for the case of D = 4 supersymmetry. 
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Figure 7: (a-b) Polarized binned angular distributions for each helicity configuration, taking 
an 80% polarization of the initial electron beam. The solid histogram represents the D = 4 
bino-like model, while the "data" points correspond to the effects of a supersymmetric bulk 
with Ac = 1.5, 3.0, 6.0 TeV from top to bottom. The Ac = 6.0 TeV case is only discernable 
in Figure (b). 
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Figure 8: (c-e) Polarized binned angular distributions for each helicity configuration, taking 
an 80% polarization of the initial electron beam. The solid histogram represents the D = 4 
bino-like model, while the "data" points correspond to the effects of a supersymmetric bulk 
with Ac = 1.5, 3.0, 6.0 TeV from top to bottom. 
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production. It is given by 



A 



pol 



dai ~ da} 
dai + da 



R 



(51) 



where the left- and right-handed subscripts refer to the polarization of the initial electron 
beam, i.e., dai = da^e^e^ e^ej , e^e^)/(icos6'. This asymmetry is displayed in Fig. ^ 
where the solid histogram again represents our D = 4 bino-like supersymmetric model and 
the "data" points are for a supersymmetric bulk with = 1.5 and 3.0 TeV. The error bars 
are again statistical only and assume an integrated luminosity of 500 fb^^. The electron 
beam polarization is taken to be 80%. We see that the asymmetry varies substantially from 
its -D = 4 value with the addition of gravitino KK exchange, thus providing an additional 
signal for a supersymmetric bulk. 
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Figure 9: Polarization asymmetry defined in the text binned in cos^. The solid histogram 
represents the D = 4 bino-like model, while the two sets of "data" points include the 
contributions from supersymmetric bulk with Ac = 1.5 and 3.0 TeV. 



We now compute the potential sensitivity to the cut-off scale from selectron pair 
production using our sample case with a bino-like lightest neutralino state for purposes of 
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demonstration. We employ the usual procedure, taking 

(-do_\ _ (_da_\ 

\ d cos 6 ' Ad \d cos 6 J 



bins 



lOd 



\ d cos t) I 



(52) 



where we include statistical errors only. We sum over both initial left- and right-handed 
electron polarization states, assuming P^- = 80%. The resulting 95% C.L. search for Ac 
from each final state, e^e^ , e^e^^j, and e^e^, is given as a function of integrated luminosity 

in Fig. for y/s = 0.5 and 1.0 TeV. We see that for 500 fb^^ of integrated luminosity, 
corresponding to design values, the search reach in the left- and right-handed selectron pair 
production channels is given roughly by Ac — 6 — 10 x ^/s, which is essentially what is 
achievable for bulk graviton KK exchange in the reaction e^e~ — * // @|- However, the 
e^e^ production channel yields an enormous search capability with a 95% C.L. sensitivity 
to Ac of order 25 x ^/s for design luminosity. This process thus has the potential to either 
discover a supersymmetric bulk, or eliminate the possibility of supersymmetric large extra 
dimensions as being relevant to the hierarchy problem. We stress that there is nothing special 
about our choice of supersymmetric parameters; our results will hold as long as selectrons 
are kinematically accessible to high energy e+e^ colliders. We conclude that selectron pair 
production provides a very powerful tool in searching for a supersymmetric bulk. 



5 Conclusions 

In summary, we have examined the phenomenological consequences of a supersymmetric 
bulk in the scenario of large extra dimensions. We assumed that supersymmetry is unbroken 
in the bulk, with gravitons and gravitinos being free to propagate throughout the higher 
dimensional space, and that the SM and MSSM gauge and matter fields are confined to a 
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Figure 10: 95% CL. search reach for in each production channel as a function of integrated 
luminosity for y/s = 0.5 and 1.0 TeV. 

3-brane. Motivated by string theory, we worked in the framework of D = 10 supergravity, 
and found that the KK reduction of the bulk gravitinos yields four Majorana spinors in four 
dimensions. We then assumed that the residual = 4 supersymmetry is broken near the 
fundamental scale M^), with only = 1 supersymmetry surviving at the electroweak scale. 

Starting with the D = 10 action for this scenario, we expanded the bulk gravitino 
into a KK tower of states, and determined the field equation obeyed by the spin-3/2 KK 
excitations. We then derived the coupling of the bulk gravitino KK states to fermions and 
their scalar partners on the brane. We applied these results to a phenomenological analysis 
by examining the effects of virtual exchange of the gravitino KK tower in superparticle pair 
production. We focused on the reaction e~^e~ e^e^ as this process is a benchmark for 
collider supersymmetry studies. Our numerical analysis was performed in the framework 
of gauge mediated supersymmetry breaking as it naturally affords a light zero-mode grav- 
itino. However, our results do not depend on the specifics of this particular model, with the 
exception of the existence of a light zero-mode gravitino state. 
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Performing the sum over the KK propagators, we found that the leading order con- 
tribution to this process arises from a dimension-6 operator, and is independent of the 
zero-mode mass. This is in stark contrast to the virtual exchange of spin-2 graviton KK 
states, which yields a dimension-8 operator at leading order. We thus found that the grav- 
itino KK contributions substantially alter the production rates and angular distributions for 
selectron pair production, and may essentially be isolated in the e^e^ channel. The resulting 
sensitivity to the cut-off scale is tremendous, being of order 20 — 25 x ^/s. 

We expect that the virtual exchange of gravitino KK states in hadronic coUisions will 
have somewhat less of an effect in squark and gluino pair production than what we have 
found here. The reason is that these processes are initiated by both quark annihilation and 
gluon fusion sub-processes, only one of which will be sensitive to tree-level gravitino exchange 
for a given production channel. The sensitivity to the cut-off scale will then depend on the 
relative weighting of the quark and gluon initial states. In addition, t-channel gravitino 
contributions will only be numerically relevant for up- and down-squark production due 
to flavor conservation; hence their effect will be diluted by the production of the other 
degenerate squark flavors and the relative weighting of the parton densities. 

Lastly, we note that virtual exchange of gravitino KK states may also have a large 
effect on selectron pair production in e^e^ collisions, which are tailor-made for t-channel 
Majorana exchanges. High energy Linear Colliders thus provide an excellent probe for the 
existence of supersymmetric large extra dimensions, and have the capability of discovering 
this possibility or eliminating it as being relevant to the hierarchy problem. 
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A Representation of the Dirac Algebra 

A convenient representation of the Dirac algebra, which simplifies the Kaluza-Klein decom- 
position, can be given in the 10 dimensional space-time as follows, 

F'^ = 7^*0 12 14, 

r3+(2,) ^ 7^ ^ as ® f5j , 

where the are standard 4 dimensional Dirac matrices, i, j = 1, 2, 3 and /i = 0, . . . ,3. Here 
a and (3 are 4x4 matrices satisfying 

{aul3j} = 0, (54a) 

and the ex's are standard pauli matrices satisfying 

aV^' = 5'^ + ie'^^a^ , (55) 

and 7^ anticommutes with the 7^^ 

{7^7'} = 0. (56) 



The a and /3 matrices can be represented as follows 



cri\ / -(T3\ /«(X2 

«i = I , "2 = , "3 = I , (57a) 

-ai / \(T3 / I zfT2. 
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A = I ° ^ I , P2 = \ ' I , Ps = \ ^ I . (57b) 

-io-2 



This representation of the Gamma matrices makes manifest the decomposition to four dimen- 
sions, since for /i = 0, . . . , 3 is the tensor product of four dimensional Gamma matrices 
with an 8 X 8 identity acting on an internal index. With these conventions, the space-like 
Gamma matrices are anti-hermitian, while the time-like Gamma matrix is hermitian. 

B Propagator for a Single Massive KK State 

To find the propagator for a single massive Kaluza-Klein state, we invert the kinetic piece 
of the operator appearing in (|33D . We solve 

for P"^, where O^'^^ is obtained by writing the Lagrangian for a single Kaluza-Klein state 
in flSSr) as 



£^ = i^^O"'^^^^. (59) 



The propagator for the mode specified by the vector of integers n is then given by 

pn 

(60) 



k"^ — mi + ie 

We are free to drop the ie convergence term as we are interested in computing t-channel 
diagrams for which t < 0, and so no poles are encountered. To solve (|5^) , we expand P^^^ as a 
linear combination of the standard set of sixteen 4x4 matrices formed from antisymmetrized 
combinations of the gamma matrices {1, 7^^, cr'^'', 7^7^, 7^}, which form a basis for complex 
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4x4 matrices, together with {rj^u, k^, ky], to generate the correct tensor structure. We then 



solve for the coefficients in this expansion. The result is 

/ h■^^l-■f \ i f k^\ f k^ \ 

P«.M. = I (^ + m„^) ^-r^^M -- 7^ + — 7^^ + — (61) 



and satisfies (on-shell) the standard condition k^P^^ = together with ''y^j^P^'^ = which 
projects out the spin | components. 

C Summation of Kaluza-Klein States 

The summation over the Kaluza-Klein states contributing to the propagator for the exchange 
of a virtual gravitino KK tower is inherently more complicated than in the case of spin-2 
exchange, and leads to a quantitatively different result. In addition, we must also include 
the effects of the finite mass of the zero- mode gravitino, mo [0. We state here the basic 
result to leading order in the cut-off Ac for the case when ^ |t|,mo, and 6 = Q extra 
dimensions, which is the example considered in the text. The result is easily generalized to 
other values of 6. 

The mass difference between neighboring KK states, \Am\ ~ l/Rc, is much smaller 
than the cut-off (|Am| ^ Ac) leading to a large number of states to be summed over. The 
mass of the individual KK state is given by ( |3T[ ) 

■ n 

(62) 



2 n ■ n 



c 



where is the common compactification radius. The near degeneracy of the KK state 
masses allows us to treat the discrete population of states on the lattice labeled by the 
vector n = [rii, . . . , ns) in the continuum limit. The number of KK states dN in the thin 
shell between m| and m| + dm'i is 

dN = p{m'l)dm\ , (63) 
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with the density function 



(47r)2r( 



2' 



The coherent sum over these states is at the amphtude level and involves terms of 
the form 



A? 



D„{t) ^ -i / dml p{ml) , (65) 



'0 



where the momentum exchange is in the t-channel with t < 0. Here, we have included an 
explicit ultraviolet cut-off Ac in order to regulate this divergent integral. As can be seen in 
Appendix 0, in the case of the gravitino propagator there are four distinct classes of Fa,H to 
consider, given by 

Fa,n = (66) 

with cr = 0, 1, 2, 3. Using (|62|) , (|66| ) and the change of variables y = ma/Xi with x = 



the integral (|65D can be put in the form 



-2z(27ri?,)V+'"' /■^^/'^ , iy^-'^)~y' 



Dait) ~ ^ / dy- ^— . (67) 

(47r)ir(f) Lo/x ' ^ + y' 

This integral can be evaluated by making use of Appell's hypergeometric function 



Fi(a; 6i, 62; c; Zi, ^2) which generalizes hypergeometric series to two variables [43|. It has the 
one dimensional integral representation 

T(a)T(c — a) 

— Fi{a]bi,b2;c;zi,Z2) = 

r(c) 

1 

m"-1(1 - uy-''-\l - uzi)-''\l - uz2y^^ du , (68) 
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cr 


iD^ (7r3i?6)-i 

<-* V C / 





A^/2 


1 


A^/3 


2 




3 


A^/5 



Table 2: Value of the integrals over the terms in the gravitino propagator in six extra 
dimensions. 



and the series expansion 



00 00 



{0')m+n{bl)m(yb2)n 

/ . 117373 ^1 ^2 ) 



m=0 n=Q 



m\n\{c)m+n 



where (a)„ is the Pochhammer symbol defined as 

{a)n = a{a + 1) . . . (a + n — 1) . 



(69) 



(70) 



For 5 = 6 extra dimensions, the result to leading order in <^ 1 is 

-t 



a 



(71) 



We note that at this order, the result is independent of the mass of the zero-mode, so long as 
this mass is much smaller than the cut-off scale. The results for the various values of cr are 
explicitly listed in Table 2. For 6 = 2 extra dimensions, the leading order result for cr = 3 is 



D 



a=3 



-2mRl Ac . 



(72) 



We note that the behavior of the sub-leading terms for 5 = 2 is quite different than that 
displayed in Table 2 for the case of 5 = 6. We do not rely on the above approximations, but 



evaluate (|67|) numerically in our analysis. 



An interesting feature of this result when applied to the gravitino propagator is that 
it modifies its qualitative structure, so that the summed propagator is dominated by a few 
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terms. This gives rise to the following {a = 3) leading order behavior (for 6 extra dimensions) 

J2p'''' - Zl^l!^ L.^ _ 1 M^^^ . (73) 



5 V 3 

n 

Following Han et al, in we take the relation between the compactification radius and 
the cut-off scale to be 

with the gravitational coupling constant k = \/WkG^. 
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